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into TTTT and KK has been performed. No clear 77c signal is observed, and upper limits for -B(r?c tttt) 
and B(r)c KK) are given at the 90% confidence level, B{J/ip jric) ■ B{rjc — ^ 7r'''7r~) < 1.1 x 10~^, 
B{J/tl) -fVc) ■ B(ric ttV) < 0.71 x 10"", B(J/^ ^ 7»?c) ■ B(?7c -» K+K-) < 0.96 x 10"^ and 
B{.J/^ -yric) ■ B{ric K'^K^s) < 0.53 x 10^^ 



Current address: Iowa State University, Ames, lA 50011- 
3160, USA 

^ Current address: Purdue University, West Lafayette, IN 
47907, USA 



Current address: Cornell University, Ithaca, NY 14853, 
USA 

Current address: Laboratoire de I'Acceleratear Lineaire, 
Orsay, F-91898, France 



2 



M. Ablikim et al.: Search for 77c Decays into tttt and KK 



1 Introduction 

Finding the source of CP violation is one of the most im- 
portant goals of particle physics. Violation of CP sym- 
metry has important consequences; it is one of the key 
ingredients for the baryon-antibaryon asymmetry in our 
universe. Early in 1964, CP violation was discovered in 
the neutral K meson system [1], and later experimental 
groups provided evidence for direct CP violation [2]. Al- 
most four decades after the original discov(Ty, CP viola- 
tion in the B meson system has boon established [3]. CP 
violation can be experimentally searched for in processes, 
such as meson decays [4] and neutrino oscillation [5], and 
by measuring the electric dipole moments of neutrons [6] , 
electrons [7] and atoms [8]. 

A total of 58 million J/ ■0 events has been collected with 
the updated Beijing Spectrometer (BESII) [9], and this 
sample offers opportunities to search for new physics in 
J/V' and r^c decays. In this ktter, a search for CP violating 
r]c decays into tttt and KK is reported using J/V' — > 77?c 
decays. 



2 The BES detector 

BESII is a conventional solenoidal magnet detector that is 

described in detail in Ref. [9]. A 12- layer vertex chamber 
(VC) surrounding the beam pipe provides track and trig- 
ger information. A forty-layer main drift chamber (MDC), 
located radially outside the VC, provides trajectory and 
energy loss (dE/dx) information for tracks over 85% of 
the total solid angle. The momentum resolution is tXp/p = 
0.017yT+p2 in GeV/c), and the dE/dx resolution 
for hadron tracks is ~ 8%. An array of 48 scintillation 
counters surrounding the MDC measures the timc-of-flight 
(TOF) of tracks with a resolution of ~ 200 ps for hadrons. 
Radially outside the TOF system is a 12 radiation length, 
lead-gas barrel shower counter (BSC). This measures the 
energies of electrons and photons over ^ 80% of the total 
solid angle with an energy resolution oioEjE = 22% /^/E 
[E in GeV). Outside of the solenoidal coil, which provides 
a 0.4 Tesla magnetic field over the tracking volume, is 
an iron flux return that is instrumented with three dou- 
ble layers of counters that identify muons of momentum 
greater than 0.5 GeV/c. 

In the analysis, a GEANT3-bascd Monte Carlo simula- 
tion program (SIMBES) with detailed consideration of de- 
tector performance is used. The consistency between data 
and Monte Carlo has been checked in many high purity 
physics channels, and the agreement is reasonable [10]. 
The detection efficiency for each decay mode is determined 
taking into account the decay angular distributions in the 
Monte Carlo simulation. The angle {6) between the direc- 
tions of e+ and r/c in the laboratory frame is generated 
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according to a 1 -|- cos^ 9 distribution, and uniform phase 
space is used for the r]c decaying into tttt and KK. 



3 Event selection 

Charged tracks are reconstructed using hits in the MDC, 
and they are required to have a good helix fit and sat- 
isfy I cos^l < 0.8, where is the polar angle of the track. 
For all decays except J/ifj ^ ^K^K^, the point of clos- 
est approach of each track to the beam line must satisfy 
^/xq+Hq < 2cm, 1 2:0 1 < 20cm, where xq and yo are the 
coordinates transverse to the beam line and zq is the dis- 
tance along the beam line from the interaction point. 

A neutral cluster is considered to be an isolated pho- 
ton candidate when the energy deposited in the BSC is 
greater than 0.05 GeV, the first hit is in the beginning six 
radiation lengths, the angle between the nearest charged 
track and the cluster is greater than 18°, and the difference 
between the angle of the cluster development direction in 
the BSC and the photon emission direction is less than 
30°. More than one photon per event is allowed because 
of the possibility of fake photons coming from interactions 
of charged tracks with the shower counter and from other 
background sources. 



3.1 TJc TTTT 

To select it^-k~ candidates, the total number of hit 

layers in the muon counter is required to be less than 
four in order to remove (7)/z events, and the sum 
of Et^+ /P-K+ and E^- /P^- is required to be less than 1 
to remove the large background from (radiative) Bhabha 
events, where E^+{E^-) and Pn.+ (P^-) are the energy 
measured in the BSC and the momentum of the 7r+(7r~), 
respectively. At least one track is required to be identi- 
fied as a pion using TOF or dE/dx information. In or- 
der to reduce background from events with 7r°, Pj^ = 
A\'Prniss\'^ sin^ 6^/2 is required to be less than 0.002 (GcV/c)^, 
where Pmiss is the missing momentum of the charged 
particles, and is the angle between the missing mo- 
mentum and the photon direction. To further suppress 
the dominant pn background, events with more than one 
photon and satisfying jm-^^ — 0.135| < 0.065 GcV /c? arc 
rejected, where m^^ is the invariant mass of two isolated 
photons. Finally, to obtain better mass resolution and to 
suppress backgrounds further, events are kinematically fit- 
ted with four constaints (4C) to the J/tp 77r"'"7r~ and 
J/ip jK^K^ hypotheses and are required to satisfy 
X^^+^- < 10 and X^^+^- < X^K+K-- there is more 
than one photon, the fit is repeated using all permutations, 
and the combination with the lowest fit is retained. The 
TT"*"??" invariant mass spectrum of events surviving selec- 
tion is shown in Fig. 1(a). No r/c signal is observed. 

For rjc tt^tt", the tt'^ mesons are identified through 
tt" 77. Four constraint kinematic fits to J/ip ^ 'jir'^Tr'^ 
are performed using all combinations of five photons, and 
the combination of photons with the smallest x^ is selected 
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and is required to satisfy X57 < 10- To select tt'^s, the 
combination with the smallest 5^a is chosen, where 

and |m-yi72 — m^o| < 0.065 GeV/c^ and Im-y^ry^ — m„o\ < 
0.065 GeV/c^ are required. To reduce background from 

events with u), events with the invariant mass of a tt*^ and 
one photon satisfying \myj^o^ ^ — 711,^1 < 0.030 GcV/cP arc 

rejected. The tt'^tt'^ mass spectrum after applying these 
selection criteria is shown in Fig. 1(b). No evident r/c is 
observed in the mass spectrum. 

3.2 KK 

The total number of hit layers in the muon counter is 
required to be less than four to remove J/ip — > jfx^ ii~ 
background. Candidate rjc K^K~ events are required 
to have P"^^ < 0.002 (GeV/c)^ to ehminate 7r° background. 
To further reduce J/ip -k^K+K' and J/tj) 
contamination, events surviving the above criteria and 
with two or more photons are kinematically fitted to these 
hypotheses, and events with a fit < 50 and with a pho- 
ton pair invariant mass within 0.065 GeV/c^ of the 7r° 
mass are rejected. Finally, to obtain better mass resolu- 
tion and to suppress backgrounds further, the two charged 
tracks and one photon in the event are 4-C kinematically 
fitted to the J/ip jn~^n~ and J/tp ^ ^K^K~ hy- 
potheses and X^^K+K- < 10 and X^^k+k- < X^tt+tt- a""® 
required. 

The K~^K~ mass spectrum is shown in Fig. 2(a). No 

rjc signal is evident. 

For the decay r?c -^^5^5, 

mesons are identified 
through —^ TT+TT . To ensure that the two charged pi- 
ous are from the same Kg vertex, \zi — z-2\ and l^^s — 2:4] 
are required to be less than 0.06 m, where z\ and 22 are 
the Z coordinates of point of closest approach of the track 
to the beam axis of the first two pions, and z^ are 
those of the second two pions. Candidate events must sat- 
isfy (5|o < (0.020GeV/c2)2, where = [m^+^-(l) - 

m^^o]^ + [i^-R+TT- (2) — Wjfo]^ and m^+jr- is calculated at 
the Kg decay vertex. The main backgrounds from jK^K^tt^ 
and jKgKgTT^ events are suppressed by requiring the 4-C 
kinematic fit x^A-n < 10- The distribution of K^Kg invari- 
ant mass is shown in Fig. 2(b). No ric signal is seen. 

4 Background estimation 

Backgrounds in the various processes except rjc KgKg 
are estimated using Monte Carlo simulation. Possible back- 
grounds and their contribution are listed in Table 1. The 
number of background event from Monte Carlo are nor- 
malized to data according to the corresponding luminosity 
(or known branching fraction) and efficiency of the sam- 
ples. The main backgrounds come from J/tp {"f)iJ,'^ii~ 
and J/tp (7)e+e^ for rjc tt+tt", from J/tp ^ uiw^ 
for ric 7r°7r°, and from J/tp ^ K*^K~ + c.c for r/c — > 



K'^K~ . The shape of the background from Monte Carlo 
simulation is consistent with that of data. 

For the decay r]c KgKg, non Kg background is esti- 
mated using the i5|.o sideband region, defined by 0.0008 < 

dl-o < 0.0012 (GeV/c2)2. The sideband histogram in Fig. 

2(b) has 18 background events in total. The shape of the 
background estimated by sideband is consistent with that 

of data. 

In fitting data, the backgrounds are fitted with poly- 
nomial function. 



5 Systematic errors 

Systematic errors on the results arise from the uncertain- 
ties in the number of J/4' events, the secondary branching 
fractions of the decay modes considered, the estimation of 
the selection efficiency, and the determination of the back- 
ground. The various contributions to the systematic error 
are listed for all decay modes in Table 2. 

1. The MDC tracking efficiency has been measured us- 
ing channels like J/ip AA and ip{2S) — » tt+tt" J/V', 
J/t/j It is found that the Monte Carlo sim- 
ulation agrees with data within 1 — 2% [10] for each 
charged track. Therefore 4% is taken as the systematic 
error for J/ijj 77r+7r~ and for J/ip — > 'yK~^K~, 8% 
for J/iP iK^K^. 

2. The systematic error on the pion identification effi- 
ciency is found by comparing the efficiency difference 
between data and Monte Carlo. It has been studied in 
detail using the decay J pn [10], where it is foimd 
that the identification efficiency for data is in good 
agreement with that of Monte Carlo. For the study of 
r]c 7r+7r~, particle identification is only required for 
one pion, and its systematic error is negligible. 

3. The systematic error from the kinematic fit has been 
estimated using J/ip — > fm [13] and J/'ip 7r"'"7r~7r"'"7r~7r' 
decays where it is possible to identify events cleanly 
without a kinematic fit. The results are 4.0% and 4.3% 
for the above two decay modes. Here 5% is taken as 
the systematic error from the kinematic fit for all decay 
modes studied in this paper. 

4. The photon detection efficiency has been studied using 
J/xp p'^Ti^ [14], and the difference between data and 
Monte Carlo simulation is about 2% for each photon. 
The resulting systematic errors on the branching frac- 
tions in this analysis range from 2% to 10% depending 
on the decay mode. 

5. The systematic error coming from fx counter simula- 
tion is estimated by studying the behavior of two pi- 
ons in the ji counter from the decays J/ip ^ p+7r~ 
and J/tp p~TT'^. The result indicates that the dif- 
ference between data and Monte Carlo simulation is 
1.2% which is taken as the systematic error. 

6. The systematic error arising from tt*^ reconstruction 
is estimated by studying J ftp p°7r°. We find the 
m-y^ distribution from data is in good agreement with 
that from Monte Carlo. The difference is about 0.4%. 
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Fig. 1. Distributions of tttt invariant mass near the rjc mass region, where (a) is for r^c —* 7r"'"7r~ and (b) is for r]c — > 7r°7r°. The 

points with error bars are data, and the histograms are Monte Carlo simulated background. The smooth curves in the plots arc 
the best fits to the data, while the shaded histograms are the expected shapes of an rjc signal as determined by Monte Carlo 
simulation (not normalized). 
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Fig. 2. Distribution of KK invariant mass near the rjc mass region, where (a) is for r?c — > K'^ K~ and (b) is for rjc — » KgKg. 
The points with error bars are data, the histogram in (a) is Monte Carlo simulated background, and the histogram in (b) is 

background estimated using the (5|,q sideband. The smooth curves in the plots are the best fits to the data, while the shaded 

s 

histograms are the expected shapes of an rjc signal as determined by Monte Carlo simulation (not normalized). 



Table 1. Dominant sources of background in the r/c tv^tt , rjc tt'^tt'^, and r/c K^ K . Hero N^ata and N^kg are the 
numbers of events and estimated backgrounds in the 2a wide rjc mass region, respectively (normalized). 



77c TT^TT 


Nbkg 




Nbkg 


rjc K+K- 


Nbkg 




585 


W7r"(727r") 


53 


K*+K- + c.c. 


7 


(7)e+e- 


281 


a;/2(1270)(737r") 


11 




5 




9 




9 


(7)/i+/i- 


5 


■yK+K- 


5 


7r/(77r'V"r7) 


5 


777 """TT" 


1 






tJ7r"7r'^(737r") 


4 






sum 


880 ± 36 




82 ±8 




18 ±1 


Ndata 


872 




84 
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Table 2. Systematic errors for the four decay modes (%). 



Sources 




u u 


r^c ^ K+K- 


77. ^ K'^sK'k 


MDC tracking 


4 




4 


8 


Photon efficiency 


2 


10 


2 


2 


Kinematic fit 


5 


5 


5 


5 


/i counter simulation 


1.2 




1.2 




tt" reconstruction 




0.9 






Mass spectrum fitting 


2.9 


4.0 


9.4 


5.6 


Background uncertainty 


11.8 


5.8 


6.5 


8.1 


Total number of J/t/j 


4.7 


4.7 


4.7 


4.7 


Total 


14.7 


14.1 


14.1 


14.6 



Table 3. Numbers used in the calculation of upper limits of branching fraction at the 90% confidence level and the upper limits 
for all modes. 



Decay modes 




e{%) 


systematic error(%) 


Br. 


J I'll) -> 777c, 7?c -> 7r+7r~ 


< 60.2 


10.9 ±0.1 


14.7 


< 1.1 X lO""" 


J/tp 7r?c,7?c 7r°7r° 


< 27.0 


7.7 ±0.1 


14.1 


< 0.71 X 10"^ 




< 30.1 


6.3 ±0.1 


14.1 


< 0.96 X 10"'' 




< 12.5 


10.1 ±0.1 


14.6 


< 0.53 X 10"-''' 


r]c tt'^tt 






34.1 


< 1.1 X 10"'^ 








33.9 


< 0.71 X 10"^ 


r?c ^ K+K- 






33.9 


< 0.96 X 10"^ 








34.1 


< 0.53 X 10"^ 



Therefore 0.9% is regarded as the systematic error for 

7. The uncertainties of the rjc's mass and width are also 
systematic errors in the upper Umit determination. 
When doing the fit to the mass spectrum, we fix the 
?7c's mass and width to their PDG values. Changing the 
r^c's mass and width by one standard deviation from 
the PDG values in fitting the mass spectrum gives sys- 
tematic errors from 2.9% to 9.4%. 

8. Another systematic error arises from the uncertainty of 
the background shape. Changing the order of the back- 
ground polynomial and the fitting range gives system- 
atic errors from 5.8% to 11.8% for the different decay 
modes. 

9. The uncertainty on the number of J/ip events intro- 
duces a systematic error of 4.7% [15], common for all 
decay modes. 

The systematic errors from the above sources are listed 
in Table 2. The total uncertainty is the sum in quadrature 
of the individual contributions. 



6 Results 

Since no evident ijc signal is observed in the four decay 
modes, we determine upper limits on the branching frac- 
tions {B) with the following formula 



B 



N. 



UL 



where N^^ is the upper limit on the number of observed 

events for the rjc signal, e is the detection efficiency ob- 
tained from Monte Carlo simulation, Nj/^ is the total 
number of J/V' events, (57.7 ± 2.7) x 10^ [15], which is 
obtained from inclusive 4-prong hadronic decays. For the 
decay ?7c K^K^, this should be divided by the square 
of the K'g TT+TT^ branching fraction, when calculating 
the branching fraction of rjc KgKg. 

Using a Breit Wigner fit to the data, the result is con- 
sistent with no signal. We determine a Bayesian 90% con- 
fidence level upper limit on Ns by finding the value N^^ 
such that 



0.90, 



-■N 



where A''^ is the number of signal events, and L is the value 
of the likelihood as a function of Ng. 

The values of A^^^, e, and the upper limits on the 

branching fractions for all decay modes are listed in Table 
3, where A^,^^ and B are given at a Bayesian 90% confi- 
dence level. The systematic errors are included by lower- 
ing the efficiencies by one standard deviation. When the 
upper limits on B(77c TT~^Tr~ ,Tr'^Tr'^ , KgKg) are 
derived in Table 3, an additional systematic uncertainty 
of 30.8% from B{J/^ 7r7c),(1.3±0.4)% [16], is included 
in the total systematic errors. 



7 Summary 

Based on 58 million J/ip events taken at BESII, a search 
for rjc rare decays into tttt and KK has been performed. No 
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clear rfc signal is observed, and upper limits for r^c tttt 
and r}c KK have been obtained at the 90% confidence 
level for the first time 

B{J/i; jrjc) ■ B{r]c ^ n+n-) < 1.1 x IQ-^, 
B{J/xP 777c) • Blr]c ^ 7r°7rO) < 0.71 x 10"^, 
B( J/V- ^ 77/0) • B{r]c K+K-) < 0.96 x 10-^ 
B(J/V' ^ 7r?c) • S(77c ^ K%K%) < 0.53 x lO'^. 
These results provide experimental limits for theoretical 
models predicting how much CP violation there may be 
in ?7c meson decays. 
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